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Abstract 

The Mel'nikov equation is a (2+1) dimensional nonlinear evolution equation ad- 
mitting boomeron type solutions. In this paper, after showing that it satisfies the 
Painleve property, we obtain exponentially localized dromion type solutions from 
the bilinearized version which have not been reported so far. We also obtain more 
general dromion type solutions with spatially varying amplitude as well as induced 
multi-dromion solutions. 



1 Introduction 



The identification of dromions which are exponentially localized solutions in 
(2+1) dimensional soliton equations [1-6] has been one of the most interest- 
ing developments in soliton theory in recent times, which has given a fillip 
to the understanding of integrable systems in (2+1) dimensions. Essentially, 
these localized solutions arise due to the presence of some additional nonlo- 
cal terms or effective local fields associated with "boundaries". Further, the 
advent of "explode decay dromions" which are again exponentially localized 
solutions with time varying amplitudes [7] and induced dromions [6,8] using 
arbitrary functions of space and time variables have set in motion the process 
of unearthing more and more novel localized entities in (2+1) dimensional 
nonlinear systems. 
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1 . 1 Mel 'nikov Equation 

An interesting evolution equation in (2+1) dimensions which we consider here 
is the one proposed by Mel'nikov [9,10] that describes (under certain condi- 
tions) the interaction of two waves on the x-axis. This equation is of the form 



3u tt - [u y + (3m 2 + u xx + 8K|x| 2 )x]a: = 0, (la) 

iXt = ux + Xxx, (lb) 

where u is the long wave amplitude (real), x is the complex short wave en- 
velope, and the parameter k satisfies the condition k 2 = 1. Eq. (1) may be 
considered either as a generalization of the K-P equation with the addition of 
a complex scalar field or as a generalization of the NLS equation with a real 
scalar field (after suitable interchange of coordinates y and t). Mel'nikov [10] 
has pointed out that Eq. (1) admits boomeron type solutions, which can be 
realized from an asymptotic analysis of the two soliton solution. 

It is expected that the investigation of this equation may have wider ram- 
ifications in plasma physics, nonlinear optics and hydrodynamics. It is this 
diverse presence of this equation which prompts one to make a detailed inves- 
tigation of their dynamics, particularly to identify whether localized solutions 
exist in this system. For this purpose, we first carry out a Painleve singularity 
structure analysis and confirm that Eq. (1) does indeed satisfy the Painleve 
property. Then bilinearizing the evolution equation and making use of cer- 
tain arbitrary functions present in the solution, we obtain a large class of 
exponentially localized dromion solutions. 



2 Singularity Structure Analysis of Mel'nikov Equation 

We explore the singularity structure of Eq. (1), by rewriting x=a and x* = D 

as 

3u tt - u xy - 6u 2 x - 6uu xx - u xxxx - 8n(a xx b + 2a x b x + ab xx ) = 0, (2a) 

ia t = ua + a xx , (2b) 
-ih t = u b + b xx . (2c) 

We now effect a local Laurent expansion in the neighbourhood of a noncharac- 
teristic singular manifold <p(x,y,t) = 0, <p x ^ 0, <p y ^ 0. Assuming the leading 
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orders of the solutions of Eq. (2) to have the form 



u = u (p a , a = a o /3 b = bo<P , 



(3) 



where Uq, ao and bo are analytic functions of (x, y, t) and a, (3, 7 are integers to 
be determined, we now substitute (3) into (2) and balance the most dominant 
terms to get 

a = p = 7 = -2, (4) 

with the condition 

Ka b = -90*, u = -6(f) 2 x . (5) 

Now, considering the generalized Laurent expansion of the solutions in the 
neighbourhood of the singular manifold 



u = u (j) a + ... + u r (j) r+a + 

< /3 + ... + a r r+/3 + .., 
b = bo^ + ... + b r (P r+ ' f + 



(6a) 
(6b) 
(6c) 



the resonances (powers) at which arbitrary functions enter into (6) can be 
determined by substituting (6) into (2) and comparing the coefficients of 
(0 r - 6 ,0 r - 4 ,0 r - 4 ) to give 



/ rl, ,:2 



V 



r — 5r — 30)f 8K,bor 
a 4>l( r2 — 5r) 

b 



Snaof 




5r) 



0, 



(7) 



where f = (r — 4)(r — 5). Solving Eq. (7), one gets the resonance values as 

r = -3,-1,0,4,5,5,6,8. (8) 

The resonance at r = -1 naturally represents the arbitrariness of the manifold 
<f)(x,y,t) = 0. In order to prove the existence of arbitrary functions at the 
other resonance values, we now substitute the full Laurent series 



u = u (/) a + Y / u r (p r+a , 

r 

r 

b = b ^ + EM> r+7 



(9a) 
(9b) 
(9c) 



into Eq. (2). Now collecting the coefficients of (</> 6 ,0 4 ,<p 4 ) and solving them, 
we obtain the relations (5), implying a resonance at r = 0. 



3 



Similarly collecting the following coefficients, we obtain the necessary infor- 
mation about the positive resonances: 



(i) coefficients of (0 5 ,0 3 ,4> 3 ): 

n ia °^* k - ib 0<t>t , nn x 
u 1 = 0,ai = — a 0x ,b 1 = — b 0x . (fO) 

(ii) coefficients of (0~ 4 ,0~ 2 ,0~ 2 ): u 2 , a 2 , and 6 2 can be uniquely determined, 
(hi) coefficients of (0" 3 ,0 _1 ,0 _1 ): u 3 , a 3 , and b 3 can be uniquely determined. 

(iv) coefficients of (0~ 2 ,</>°,0°): Only two equations result for three unknowns 
«4, a 4 , 6 4 and so one of them is arbitrary, corresponding to a resonance at r 
= 4. 

(v) coefficients of (4>~ l ,<f> 1 ,4> l ): Only « 5 is determined, while a 5 and 65 are 
arbitrary corresponding to double resonance at r = (5, 5). 

(vi) coefficients of (0°,0 2 ,0 2 ): Only two equations result for three unknowns 
u 6 , a 6 , b 6 and so one of them is arbitrary, corresponding to r = 6. 

(vii) coefficients of (0 1 ,0 3 ,0 3 ): u 7 , aj, and b 7 can be determined in terms of 
earlier coefficients. 

(viii) coefficients of (c/> 2 ,(/> 4 ,0 4 ): Only two equations result for three unknowns 
u 8 , a 8 , b 8 and so one of them is arbitrary, corresponding to r = 8. 

For the negative resonance r = -3, following the approach of Conte, Fordy and 
Pickering [11], we demand that both the solution of Mel'nikov Eq. (2) and the 
solution close to it represented by a peturbation series in a small parameter 
e are free from movable critical manifolds. We identify that the first order 
peturbed series does admit an arbitrary function corresponding to a movable 
pole at the resonance r = -3. Consequently, for each of the eight resonances 
given by Eq. (8), one can associate an arbitrary function in the solution (and 
close to it) without the introduction of movable critical manifolds. 

It must be mentioned that the above system (2) admits another leading order 
behaviour with a = — 2, [3 = 7 = — 1, uq = — 20 2 , ao and bo are arbi- 
trary. The Laurent series with the above leading order leads to resonances 
0,0,-1,3,3,4,5,6, corresponding to a normal branch and the existence of 
sufficient number of arbitrary functions can be established at these resonance 
values. This can also be verified from the corresponding bilinear form as stud- 
ied by Grammaticos, Ramani and Hietarinta [12]. Consequently one can be 
assured that the Mel'nikov equation (1) or (2) indeed satisfies the Painleve 
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property. 



3 Bilinearization of Mel'nikov Equation and Localized Solutions 



We next Hitora bilinearize the Mel'nikov Eq. (1) to bring out the existence of 
exponentially localized solutions of Mel'nikov equation. Making the transfor- 
mation 



q2 

u = 2— InF, (11a) 
X = % (Hb) 

identifiable from the Painleve analysis, Eq. (1) gets converted into the follow- 
ing Hirota bilinear form, 



(3D 2 - D x D y - Dl)F.F = 8k\G\ 2 , (12a) 
iD t G.F = D 2 X G.F, (12b) 

where the D's are the usual bilinear operators. However, this bilinearization 
have been already done by Y. Hase etal [13] where they have given soliton solu- 
tions whereas we have brought out localised solutions here. For a completion, 
we proceed as follows. Introducing the series expansion, 



G = egW + (?gW + ..., (13a) 

F = l + e 2 / (2) +e 4 / (4) + _ (13b) 

into the above bilinear form and gathering terms with various powers of the 
small parameter e, we obtain the following set of equations, 



0(e):id l) =gW, (14a) 
0(6 2 ):3/f - - /Ex = WW (14b) 

etc. Solving Eq. (14a), we can immediately write down the following solution, 



= ^2 exp(ipj),ipj = IjX + rrijy + Ujt + ipf\iuj = I 2 , (15) 
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where the spectral parameters lj,rrij,ujj and tp^ are all complex. Confining 
to N = 1 in Eq. (15) and substituting (15) into (14b), we obtain 

/W = ^ + * + 2 ' 4 >- CXp < 2 - 4 » = (3^-^-4^) - < 16 > 

Here l\ = l\ R + ilu, mi = mm + imu and uj\ = u>i R + iuju and also iu)\ = l\. 
Choosing g&+V = 0, /(«) = 0, for j > 1, in Eq. (13) and using Eqs. (15) and 
(16) alongwith the transformation (11), the physical field u and the potential 
X can be easily seen to be driven by the envelope soliton and pulse soliton 
respectively as 



x = ( ^(i^^-^n sechW , is + (17a) 

u = 2l 2 1R sech 2 (tp 1R + A). (17b) 

where ip 1R = l 1R x + m 1R y + 2l 1R l u t + ipfl and ip u = l u x + m u y - (lf R - 
can proceed further in the standard way to obtain higher 
order soliton solutions also. 



3. 1 Dromions 

Looking at the above solutions (17), we realize the fact that as the parameter 
^lr 0, both u and x vanish. But, when ni\ R — > VH\ R l\j — 4lf R , the potential 
X vanishes, whereas the physical field u survives and is driven by a ghost 
soliton of the form, 

M = 2/^sech 2 (^ li? + c), (18) 

where c is a new constant. This predicts the existence of exponentially localized 
solution for the complex field variable x 1X1 Eq. (1). 

3.1.1 (1,1) Dromion 

To generate a (1,1) dromion, we now make the ansatz 

p — i _|_ gV'iH-^ _|_ e i>2+i>2 _|_ K e i>i+i>t+^2+^l ^ ^) 

where 

vf)i=px + ut, (20a) 
i>2 = qy- (20b) 
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Here K is a real constant and p, u and q are complex constants. Substituting 
(19) in (12a), we obtain 



G = pe^ 1+ ^ 2 , (21a) 
4k\ P \ 2 = (p + p*)(q + q*)(l - K) (21b) 

for the parametric choice 3u; R = Ap R . Substituting (19) and (21) in (12b), we 
find that u> = —ip 2 , where p R = ±^/3pj, pr and pj are the real and imaginary 
parts, respectively of p. Hence, the exponentially localized solution with one 
bound state for the potential field x f° r the above choice of u> and p takes the 
form 

X= — — - - - , (22) 

while the scalar field has the form 

4p2(i + e ^+^)( e ^i+^ + # e &W+&+#) 

w = = — = = — = = — = — = — = , (23) 

(1 _|_ e V>i+^ _|_ e V-2+V| _|_ Ke^+W+^+W) 2 

which is always bounded, but localized everywhere except in the neighbour- 
hood of the line ip± + = in the (x,y) plane. A snapshot of the (1,1) 




dromion solution for the magnitude of the potential field \ is shown in Fig. 1. 
One can proceed to find multi-dromion solutions also, generalizing the above 
(1,1) dromions. 
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3.1.2 Dromions with Spatially Varying Amplitude 

It can be seen from Eq. (14a) that the differential equation for g involves only 
the variables t and x. Hence, an arbitrary function of y can also enter into its 
solution so that the most general form of it can be given as 

N 

g W = exptyj),^ = l jX + fj(y) + Wjt + ^■\iw j = l], (24) 
j=i 

where fj(y)s are arbitrary functions of y. This fact can be harnessed in a 
suitable way to construct a more general class of localized solutions. Following 
the above procedure to derive the (1,1) dromion solution (22), one can easily 
obtain the generalized dromion solution involving arbitrary function of y in 
the same form as Eq. (22) except that ip 2 is now given by 

h = f(y), (25) 

instead of (20b), where the arbitrary function f(y) is in general complex. The 
amplitude of the localized solution with arbitrary function of y is now defined 
by the equation 

AK\p\ 2 = 2(p + p*)(l-K)f' R (y). (26) 

where f' R (y) is the derivative of the real part of / with respect to y. Thus, the 
amplitude of the above localized solution varies with the spatial coordinate 
y by virtue of Eq. (26). This situation is reminiscent of the explode-decay 
dromion of the variable coefficient DSI equation [7] where the amplitude varies 
with time. It should be mentioned that, to our knowledge this is the first 
time the amplitude of a localized solution of a (2+1) dimensional nonlinear 
partial differential equation has been found to vary as a function of the spatial 
coordinate y. This can be easily generalized to multi-dromions with spatially 
varying amplitude. 



3.1.3 Induced dromions 

We also wish to point out that the existence of an arbitrary function in the 
solution of g^ in Eq. (14a) can be further utilized to obtain new induced 
dromion solutions [6] for \- F° r example, Eqs. (14a) and (14b) can also be 
solved in terms of arbitrary functions as 

g W =a(y)e^J 1 = l 1 x + u; 1 t, (27) 

where a(y) is an arbitrary complex function of y and l\ and u)\ are complex 
constants constrained by the condition %lo\ = l\. Substituting the form (27) 
into (14b) and solving the resultant equation, we obtain 

/( 2 ) = b(y)e^ + ^*, (28) 
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where b(y) is a real function of y and is given by the condition 
Qu 2 1R b(y) - h R b y - 8l 4 1R b(y) = 2K\a(y)\ 2 . 



(29) 



The above solutions can then be used to generate curved line soliton for the 
field variable and the potential as 



X = 



2^ 



e"^secl# 1R + ^log%)], 



u = 2l 2 lR sech 2 [if 1R + -log%)]. 



(30a) 
(30b) 



where ip\ R and ipu are the real and imaginary parts, repectively of ipi. Thus, by 
choosing the arbitrary functions a(y) and b(y) suitably which are constrained 
by the Eq. (29), one can induce localized solutions for the field variable x as 
in the case of Zakharov-Strachan equation [6]. Eventhough there exists two 
functions a{y) and b(y), only one of them is found to be arbitrary, which is 
evident from the Eq. (29) above. For example, by choosing 



g(y) 



2sech(mi2/), 



where m\ is a real constant, we can find from Eq. (29) that 

b(y) = exp f--^-(8«;tanhmiy + (8/^ - 6w 2 R )y) 

Then we obtain the induced localized solution 

1 



(31) 



(32) 



X = e* v ' 1/ sech(m 1 y)sech 



l 1R (x + 2l u t)- 



9/ r ^8/itanhmiy + (8l 1R - QuJ 1R )y 



(33) 



Thus, by choosing a(y) and b(y) suitably, one can induce a wide class of local- 
ized solutions for the Mel'nikov Eq. (1). For example, choosing an algebraic 
form 

= I (34) 

we obtain an algebraically decaying localized solution 



1 



X = 



(y + y ) 2 + 1 



^"sech 



l 1R (x + 2lut) 



21 



1R 



(Sl 4 1R - 6u 2 1R )y 



+8k 



I 



dy 



((y + y ) 2 + l) 



(35) 
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One can as well generalize this procedure to construct even wider class of 
localized solutions. In fact, multi-induced dromions take the simple form 



Xn 



2^ 



e^sech[^ 1R + -log%)], 



(36) 



where a,j, j = 1, 2...7V are arbitrary functions of y and they are related to b(y) 
by the relation 



N N 



^IrKv) - hnby - 8l 4 1R b(y) = 2kJ2J2 a p{v) a * q (y)- 



(37) 



p=l q=l 



Note that the sum of the arbitrary functions of y on the right hand side of 
eq.(36) becomes possible due to the structure of eqs. (1). In practice one can 

choose ^ J= ^— — ^ - conveniently, for example, as a combination of algebraic and 



(38) 



hyperbolic functions. With N=2, choosing the functions as 

! ' - 2sech(miy + 81), a2 ^}_ = 2sech(m 2 y + 5 2 ), 



where mi, m 2 , Si and 82 are parameters, solving eq.(29), one can obtain 



Ky) = exp - —((8lf R - Qujl R )y + 8/<tanh(mii/ + 81) 
\ i>ir 

+tanh(m 2 y + 8 2 ) + 2 J sech(miy + 8i)sech(m 2 y + 8 2 )dy)) \ . (39) 



Then, the induced two dromion solution is given by 



X2 = e^ 17 (sech(mi2/ + Si) + sech(m 2 y + 5 2 ))sech l 1R (x + 2lut) 



(8li R — Quo\ R )y + 8«(tanh(miy + Si) 



21 



1R 



+tanh(m 2 y + 5 2 ) + 2 J sech(miy + <5i)sech(m 2 j/ + 8 2 )dy) 



(40) 



and is shown in Fig. 2. One can identify the mutual influence of one dromion 
over the other from the additional terms occurring in the square bracket of 
eq. (40). 
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Fig. 2. Snapshot of the induced two dromion solution of the Mel'nikov equation (see 
Eq. (40)) ' 



If one chooses both the functions cij, j = 1,2, algebraically, 
ai(y) 2 a 2 (y) 2 



y / %) {m 1 y + y w ) 2 + 1' y Jb(yj ( m ^V + Vw) 2 + 1 
or one of the functions to be algebraic and the other one to be hypebolic 



(41) 



= 2sech(m l|/ + 5,),^ = - ( 2 —— (42) 

one can generate different kinds of induced lump-lump or dromion-lump solu- 
tions, respectively. 

We have also tried to generate more general two soliton solutions by choosing 

g {1) = ai^)^ 1 +a 2 {y)e^\i) l = hx + 10^,^2 = hx + utf, (43) 

which lead to a condition li = l 2 thereby reducing to our original form (36) 
for N = 2. This is also true for N > 2. Thus we believe that the solution 
(36) constitutes the most general localized solution we could construct for the 
Mel'nikov equation through our procedure. 



4 Conclusion 



In this paper, we have pointed out the interesting fact that the (2+1) dimen- 
sional Mel'nikov equation admits exponentially localized solutions of different 
classes. We have also checked its integrability through Painleve analysis. We 
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have in particular constructed localized dromion solutions and obtained new 
classes of localized solutions such as dromions with spatially varying amplitude 
and induced dromions. 
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